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We present molecular dynamics simulations of a simple model for polymer melts with intramolecu- 
lar barriers. We investigate structural relaxation as a function of the barrier strength. Dynamic cor- 
QQ , relators can be consistently analyzed within the framework of the Mode Coupling Theory (MCT) of 

■ the glass transition. Control parameters are tuned in order to induce a competition between general 

' packing effects and polymer-specific intramolecular barriers as mechanisms for dynamic arrest. This 

^Nj , competition yields unusually large values of the so-called MCT exponent parameter and rationalize 

qualitatively different observations for simple bead-spring and realistic polymers. The systematic 
O ' study of the effect of intramolecular barriers presented here also establishes a fundamental difference 

I between the nature of the glass transition in polymers and in simple glass-formers. 

CN . PACS numbers; 64.70.pj, 64.70.qj, 61.20.Ja 
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Since they do not easily crystallize, polymers are 
probably the most extensively studied systems in re- 
lation with the glass transition phenomenon. Hav- 
ing said this, their macromolecular character, and in 
particular chain connectivity, must not be forgotten. 
Its most evident effect is the sublinear increase of the 
mean squared displacement (Rouse-like) arising after 
the decaging process, in contrast to the linear regime 
found in non-polymeric glass-formers. Another particu- 
lar ingredient of polymers is that, apart from fast libra- 
tions or methyl group rotations, every motion, as local 
as it be, involves jumps over carbon-carbon rotational 
barriers and/or chain conformational changes. 

In this work we investigate, by means of molecular 
dynamics simulations, the decisive role of intramolec- 
ular barriers on the glass transition of polymer melts, 
by systematically tuning barrier strength in a simple 
bead-spring model. We discuss the obtained results 
within the framework of the Mode Coupling Theory 
(MCT) of the glass transition [H. Initially derived for 
monoatomic hard-sphere systems, the theory has been 
further developed for more complex systems, includ- 
ing fully-flexible bead-spring chains as simple models 
for polymer melts . MCT asymptotic laws have been 
tested in different polymeric systems. The values of the 
associated dynamic exponents exhibit significant dif- 
ferences between the limits of fully-flexible bead-spring 
chains 3] and fully-atomistic polymers ^^]. In partic- 
ular, the so-called exponent parameter takes standard 
values A ~ 0.7 for the former case and values approach- 
ing the upper limit A = 1 for chemically realistic poly- 
mers [3] . While the former A- values are characteristic of 
systems dominated by packing effects, as the archetype 
hard-sphere fluid, the limit A = 1 arises at higher-order 
MCT transitions The latter, or more generally tran- 



sitions with A ^ 1, arise in systems with different com- 
peting mechanisms for dynamic arrest. These systems 
include short-ranged attractive colloids 0, [3| (competi- 
tion between short-range attraction and hard-sphere re- 
pulsion) or binary mixtures with strong dynamic asym- 
metry [^j Q (bulk- like caging and confinement). 

Motivated by these analogies, we argue that values 
A ^ 1 for real polymers also arise from the compe- 
tition between two distinct mechanisms for dynamic 
arrest: usual packing effects and polymer-specific in- 
tramolecular barriers. Such barriers are not present 
in fuUy-fiexible bead-spring chains, which exhibit stan- 
dard A- values In order to shed light on this ques- 
tion, we perform a systematic investigation of the in- 
terplay between packing and intramolecular barriers. 
Starting from fully-flexible bead-spring chains, stiffness 
is introduced by implementing intramolecular bending 
and torsion terms. The barrier strength is systemati- 
cally tuned in order to induce competition between the 
former two mechanisms. We restrict to stiffness for 
which no orientational order is present, and provide a 
complete dynamic picture of the isotropic phase as a 
function of the barrier strength. An extensive test of 
MCT asymptotic laws is performed. Simulation results 
are described with consistent sets of MCT exponents. 
A progressive increase of A is induced by strengthening 
the competition between packing and intramolecular 
barriers, confirming the proposed scenario. 

We simulate bead-spring chains of iV = 10 identical 
monomers of mass m = 1. Monomer- monomer inter- 
actions are given by a corrected soft-sphere potential: 
V{r) = 4e[(CT/r)i2 - Co + C2{r/a)^], where e = 1 and 
(7=1. V{r) is set to zero for r > ccr, with c = 1.15. 
The values Co = 7c^^^ and C2 — 6c~^^ guarantee 
continuity of potential and forces at r = ca. V{r) 
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is purely repulsive and has no local minima. Thus, 
it drives dynamic arrest only through packing effects. 
Chain connectivity is introduced through a FENE po- 
tential i, VpEMr) = -eKpRlHl - {Roa)-^r% 
between consecutive monomers, with = 15 and 
Rq = 1.5. We implement intramolecular barriers by 
means of the bending, Vb, and torsion potential, Vt, 
proposed in Ref. [l3| (see discussion there), which are 
defined for each i-momomer (1 < i < N) as: Vsidi) = 
{eKB/2){cos6i - cos6'o)^, and VT(6'i, 6*^+1, (/'i.i+i) = 
eKt s'm^ 6i s\v? 6i+i J2n=o 4>i,i+i ■ Chain stiff- 

ness is tuned by varying Kb and K^. Oq — 109.5° and 
9i is the bending angle (for 2 < i < A^ — 1) between con- 
secutive monomers i — 1, i, and i + l. The dihedral angle 
4>i,i+i is defined for the consecutive monomers i — 1, i, 
i + 1, and i + 2 (2 < i < iV — 2), as the angle between 
the two planes defined by the sets {i — 1, i, i + 1) and 
(i, i + 1, i + 2). The values of the coefficients a„ are 
0: ao = 3.00, ai = -5.90, Oa = 2.06, and = 10.95. 

Temperature T, time t, wave vector q, and monomer 
density p are given respectively in units of e/fcs (with 
ks the Boltzmann constant), (T(m/e)^/^, cr~^, and cr"'^. 
We investigate, at fixed p = 1, the T-dependence of the 
dynamics for different values of the bending and torsion 
strength, (Kb^Kt) = (0,0), (15,0.5), (25,1), (25,4), and 
(35,4). The case {Kb,Kt) = (35,4) is also studied for 
p — 0.93. The total number of chains is Nc — 300. Peri- 
odic boundary conditions are implemented. Equations 
of motion are integrated in the velocity Verlet scheme 
[m . The system is prepared by placing the chains ran- 
domly in the simulation box, with a constraint avoiding 
monomer core overlap. The initial monomer density is 
p = 0.375. Equilibration consists of a first run where 
the box is rescaled periodically by a factor 0.99 < / < 1 
until the target density p is reached, and a second iso- 
choric run at that p. Thermalization at the target T is 
achieved by periodic velocity rescaling. Once the sys- 
tem is equilibrated, a microcanonical run is performed 
for production of configurations, from which observ- 
ables are computed. For each state point, the latter 
are averaged over typically 40 independent samples. 

Orientational ordering (induced by chain stiffness) is 
discarded for all the analyzed cases by measuring the 
quantity P2(6) = (3(cos2 9) - l)/2, where 9 is the 
angle between the end-to-end vectors of two chains, and 
average is performed over all pairs of distinct chains. In 
all cases we obtain negligible values |P2(9)| < 10^^. 

We compute density-density correlators, defined 
as F{q,t) = {p{q,t)p{-q,0))/{p{q,0)pi-q,0)) where 
p{q,t) = exp[iq • rj{t)], the sum extending over 
the positions rj of all the monomers in the sys- 
tem. Density self-correlators are defined as Fs{q,t) — 
(iViVc)-iSj exp{iq- [r^ (i) - (0)]}. Results for the for- 
mer quantities are shown in Fig. [U at several g- values, 
for two state points with non-zero barriers, at T close to 



the critical MCT temperature (see below). As usual, a 
plateau is observed in the interval corresponding to the 
caging regime, i.e., the temporary trapping of a particle 
by its neighbors. This interval is known as the /3-regime 
within the framework of MCT. The second decay, cor- 
responding to full relaxation of density fluctuations of 
wave vector q, is known as the a-regime, and is often 
described by an empirical Kohlrausch- Williams- Watts 
(KWW) function, A, exp[-(t/r^K)''«], where Ag, the 
KWW time and the exponent jSq are q-dependent. 

Next we summarize the basic predictions of MCT and 
test them in the present system. In its ideal version, 
MCT predicts a sharp transition [l[ from an ergodic hq- 
uid to an arrested state (glass) at a given value of the 
relevant control parameters — here x = (T, p, Kb, Kt). 
When crossing the transition point x = Xc the long- 
time limit of F(q,t) and Fs{q,t) jumps from zero to a 
non-zero value, denoted as the critical non-ergodicity 
parameter (/^ and Z^", respectively). MCT predics 
asymptotic laws for dynamic observables. Such laws 
are characterized by dynamic exponents that are q- and 
state-point independent. They are univoquely deter- 
mined by the static correlations at x = Xc [H- More- 
over, all the dynamic exponents are univoquely related 
to a single one, the exponent parameter A (see below), 
which is the basic one controlling all MCT asymptotic 
laws. Now we summarize the main ones. 

For ergodic states close to Xc, the initial part of the 
a-process (i.e., the von Schweidler regime) is given by 
a power law expansion [ij : 



F{q,t)^rg~hg{t/T^)'> + h^''\t/T^) 
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(and analogously for self-correlators) with < 6 < 
1. The non-ergodicity parameters and the prefactors 

(2) 

hq and hq only depend on q and are different for 
each correlator. The a-relaxation time Tq, only de- 
pends on the separation parameter |x — Xc|. MCT 
predicts a divergence 12| according to the power law 
Tq oc |x — XcI"'''. In practice can be defined as 
the time where i^((j'max, decays to some small 
value z far below the plateau, with gmax the q-value 
at the maximum of the static structure factor S{q) — 
{NNc)-^{p{q,0)p{-q,0)}. Here we will use To.2- The 
exponent 7 is given by [l| : 



7=(l/2a) + (l/26). 



(2) 



As mentioned above, the full a-decay can be described 
by a KWW function. In the limit of large q MCT pre- 
dicts [13] for the KWW times a power law cx g"^/**. 
The exponents a, b, and 7 are univoquely related to the 
exponent parameter A < 1 through [l|: 



A = 



r^(i + 6) _ r^(i-a) 

r(l + 26) ^ r(l - 2a)' 



(3) 
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FIG. 1: (color online) Symbols: simulations results, at p = 
1, for density correlators. Panel (a): F{q,t) for Kb = 15, 
Kt = 0.5, at T = 0.81. Panel (b): F^{q,t) for Kb = 35, 
Kt = 4, at T = 1.33. Identical symbols in both panels 
correspond to identical wave vectors q [values are given in 
panel (a)]. Lines are fits to Eq. ([TJ. 
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FIG. 2: (color online) Critical non-ergodicity parameters, 
as determined from fits to Eq. ([l|, for different barrier 
strength at p = 1. Lines are guides for the eyes. 



with r the Euler's Gamma function. 

When numerical solutions of the MCT equations are 
not available, the former non-ergodicity parameters, 
prefactors and exponents are obtained as fit parame- 
ters from simulation or experimental data. Consistency 
of the analysis requires that dynamic correlators and 
relaxation times are described by a common set of ex- 
ponents, univoquely related through Eqs. I|2l3p . This 
consistent test has been done for all the systems here 
investigated, with different strength of the intramolec- 
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FIG. 3: (color online) Symbols; for p = 1, T- and q- 
dependence of relaxation times (see text for notations) for 
different barrier strength. Lines in (b) and (c) are fits to 
respectively (x q~^^^ and cx (T — Tc)"'' (see text). 



ular barriers. Figs. [T][3] display some representative 
examples. Fig. [1] shows at fixed p = 1 and for a broad 
g-range, fits to Eq. ([T]) of density correlators for the 
state points Kb = 15, Kt = 0.5, T = 0.81 (SI) and 
Kb = 35, K'T^A,T= 1.33 (S2). A good description 
is achieved, for all the g-values and over several time 
decades, with a fixed 5-exponent (6 = 0.50 and 0.37 
for respectively SI and S2). Fig. [2] displays, for the 
former barrier strength, the g-dependence of the crit- 
ical non-ergodicity parameters. The fully flexible case 
Kb = Kt = is also included. As deduced from the 
stronger decay of and f^^ for stronger barriers, chain 
stiffness induces a weaker localization at fixed density. 
By making an approximate fit of f^^ to Gaussian be- 
havior, exp(— 5^/^/6), we estimate, at fixed p = 1, a 
localization length Ic = 0.19, 0.21, and 0.23 for respec- 
tively {Kb,Kt) = (0,0), (15,0.5), and (34,4). 

The increase of the barrier strength at fixed p also 
induces a higher critical temperature Tc, and a longer 
relaxation time for fixed p and T. This is demonstrated 
in Fig. [31 which also shows a test of the predictions 
T„ (X (T - Tc)-'^ and oc q'^/'' (for large q) for the 
former values of the barrier strength. A good descrip- 
tion is obtained with the same ^-exponents used for the 
von Schweidler fits of Fig. [1] and with the 7-exponents 
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TABLE I: Values of the MCT exponents and critical tem- 
perature Tc for different p and barrier strength. Also in- 
cluded are the mean chain end-to-end radius Ri^ at Tc. 



System 



7 A 



Hard spheres 
Orthoterphenyl 
Polyethylene (UA) 
1,4-Polybutadiene (UA) 
1,4-Polybutadiene (FA) 
Poly (vinyl ethylene) (FA) 0.18 0.24 4.9 0.93 



0.31 0.58 2.5 0.74 

0.30 0.54 2.6 0.76 

0.27 0.46 2.9 0.81 

0.21 0.30 4.1 0.90 

0.18 0.24 4.9 0.93 



transitions. The systematic study presented in this 
work strongly suggests a competition between gen- 
eral packing effects and polymer-specific intramolecular 
barriers as the origin of this difference, ft also suggests 
a fundamental difference in the nature of the glass tran- 
sition in real polymers — driven by the former compet- 
ing mechanisms — as compared to simple glass-formers 
[ij]. Real polymers are thus classified in the family 
of complex systems as short-ranged attractive colloids 
0] or binary mixtures with strong dynamic asym- 
metry 0, 0] , which are characterized by an underlying 
higher-order MCT transition — or at least by unusually 
large values of A — arising from a competition between 
distinct mechanisms for dynamic arrest. Finally, results 
reported here provide fundamental information for mi- 
croscopic theories (and in particular for MCT) of the 
glass transition in polymers, which need to account for 
the decisive role of intramolecular barriers. 

We acknowledge support from NMP3-CT-2004- 
502235 (SoftComp, EU), MAT2007-63681 (Spain), 
2007-601021 (Spain), and IT-436-07 (GV, Spain). 



TABLE II: MCT exponents for different glass-formers. 
Data are taken from [J| and references therein. UA and FA 
denote respectively coarse-grained united atom and fully- 
atomistic models. 

derived from them through Eqs. (|2l3p . This demon- 
strates the consistency of the data analysis. For com- 
parison. Fig. [3] also includes results for the fully flexible 
case Kb = Kt = 0. 

Table H] displays the results of the MCT analysis (dy- 
namic exponents and Tc) for all the investigated cases. 
It also includes the mean chain end-to-end radius at Tc, 
i?ccj computed from the simulations. i?cc provides a 
qualitative characterization of chain stiffness. From nu- 
merical values in Table[T]a clear correlation between the 
exponent parameter A and chain stiffness is unambigu- 
ously demonstrated. The competition between packing 
effects and intramolecular barriers induces a progressive 
increase of A from the value A = 0.76 for fully- flexible 
chains to A = 0.89 for the stiffest investigated chains. 

This observation rationalizes the large difference ob- 
served between MCT exponents for fully-flexible bead- 
spring chains and chemically realistic polymers. Table 
ITTl shows a representative compilation of exponents for 
glass-formers of very different nature. Exponents for 
fully-flexible bead-spring chains are similar to those of 
non-polymeric glass-formers, including the hard-sphere 
fluid, i.e., the archetype glass- former dominated by 
packing effects. Chemically realistic polymers of in- 
creasing complexity exhibit instead values approaching 
the limit A = 1 characteristic of higher-order MCT 
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